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Abstract. We present the current performances of the AMBER / VLTI instrument in terms 
of differential observables (differential phase and differential visibility) and show that we are 
already able to reach a sufficient precision for very low mass companions spectroscopy and mass 
characterization. We perform some extrapolations with the knowledge of the current limitations 
of the instrument facility. 

We show that with the current setup of the AMBER instrument, we can already reach 3cr = 
10 -3 radians and have the potential to some low mass companions characterization (Brown 
dwarves or hypothetical very hot Extra Solar Giant Planets). With some upgrades of the VLTI 
infrastructure, improvements of the instrument calibration and improvements of the observing 
strategy, we will be able to reach 3cr = 10~ 4 radians and will have the potential to perform 
Extra Solar Giant Planets spectroscopy and mass characterization. 

Keywords, techniques: interferometric, stars: planetary systems, instrumentation: high angular 
resolution 


1. Introduction 


In this paper we discuss the current highest performances of Colour-Differential Inter¬ 
ferometry (CDI) on the AMBER instrument and compare these performances to signal 
amplitude we computed from low-mass companion simulations. 

CDI is based on simultaneous interferometric observations in different spectral chan¬ 
nels. As a high-angular resolution and high-dynamic technique, it presents two major 
advantages. First, the chromatic differences in visibility and phases are much less sensi¬ 
tive to instrumental and atmospheric instabilities, and therefore are easier to calibrate 
than the absolute complex visibility. Since the beginning of long-baseline optical inter¬ 
ferometry with separated aperture s, many early astrophysical results have been obtained 
using this self-calibration feature I Thom et al .11 198(1 iMourard et inmij. Second, the 


colour-differential phase can be measured with an accuracy much better than the angular 
interferometric resolution A /B. For objects much smaller than the diffraction limit, it is 
proportional to the variation of the object photocentre with wavelength. 

This paper is placed in the context of Extra Solar Planets characterization with inter- 
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ferometry l Yannier et al and is intending to show that this technique has already 

the potential to get some scientific results on high contrast binaries. 


2. Differential Observables computation and error bars 


2.1. The differential phase estimator 
The expression of the coherent flux on tha AMBER instrument is 


I Millour et ablEofM - ): 


N x 


C(t, A) = 2N(X)Vi(t, X)V(X)y/ Pl (t,X) P2 (t,X) 


N 


^ a ifc( A )«2fc(A)xe^ (t ’ A)+ ^ (t ’ A)+ ^ (t ’ A)+ ^ (t ’ A) 
k =1 


( 2 . 1 ) 

In this equation k is the pixel index (spatial direction), N( A) is the unknown object’s 
flux, P i(t, A) and p 2 (t, A) are transmission coefficients for the two combined beams, and 
aifc(A) and a 2 fc(A) are related to specific features of each pixel (shape of the beam). 
V)(t, A) is the instrumental contrast and V(X) is the amplitude of the complex visibility. 

(f>p(t, A) is the phase induced by the achromatic piston, (f>i(t, A) is the instrumental- 
induced phase that varies with time (since the fixed part is alread y removed by the 
data reduction algorithm of AMBER as explained in iMillonr et a,l.l ( 2004 1'). cf> 0 (t, A) is 
the observed object’s phase and <j> c (t, A) is the chromatic phase induced by other causes 
(Chromatic atmospheric phase for example). 

If we suppose that A) = 0 (no variable instrumental phase), <j> c (t,X) = 0 (no 
chromatic effect) and (f> 0 (t,X) = 0 (unresolved or centro-symmetric object) then we can 
correct the complex coherent flux from the achromatic piston effect by: 


Cnop(f,A) =C(t,X) x (2.2) 

One can note that we need an estimation of the achro matic p iston for each sample 
of time. This is a part of the Ph.D. thesis of Eric Tatulli i TatnllfEiol and it will not 
be explained in detail in this article. We compute a reference channel for each spectral 
channel, taking care of non biasing the interspectral term by removing the “work” spectral 
channel before averaging: 


A*,) = (C nop (t, Aj))^^ (2.3) 

We compute then the interspectral term between this reference spectral channel and 
the work spectral channel: 


W( A fc ) 


Cnop(t, Afc)C re f(f, Afc) 4 

ICW(t, Afc)| 2 


Then we compute the differential phase: 


(*) 


(2.4) 


^diff(A) = arg(W(A)) (2.5) 

This expression of the differential phase is quite accurate with a good achromatic piston 
correction. That is why it can be applied only to high flux sources (for example the star 
51 Peg. has a K magnitude of 5, which is sufficient for this type of application). We can 
express the interspectral term by: 


W{ A) = |W(A)|e i<?idiff(A) 


(2.6) 
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Please note that when <f> 0 (t,\) cannot be neglected with regards to A), then 
the expression of the differential phase is different and contains a bias related to the 
interferometric phase, which leads to a bias in W{ A). For well resolved objects, this effect 
has to be taken into account and leads to a specific treatment. Thus it is well beyond 
the scope of this paper and will be discussed in a futher one. 

2.2. The differential visibility estimator 

Going back to the interspectral term evaluation 12.41 we see that its modulus can be 
expressed as: y* ■ ^ we have an unbiased estimate of this modulus (which we call 

Vdiff(A)), then we can compute the differential visibility. This estimate can be made with 
the real part of the differential phase correction of the interspectral term: 

V^g(A) = 3t(W(A) x e -^ d ‘«(A)) ( 2 . 7 ) 

This estimate of the differential visibility is unbiased as we expect the visibility and 
phase to be uncorrelated. 


2.3. The error bars 

2.3.1. In theory 

Starting from the theoretical estimation of the differential phase, we can express the 
differential phase and visibility nois es f rom the f undamental photon VTV*, thermal y/N t h 
and detector orc>n y'ni n p j x noises ( Petrov 1989 1. 


(7$ = 


\J(N * “h Nth “t” ^pix ^RON ) /^ 

vW) 


a v = 


\/ N* + TVth + nf n p ix cr 


RON 


(N) 


For information, the closure phase noise is given by: 


( 2 . 8 ) 

(2.9) 




= V3 


X (7$ 


( 2 . 10 ) 


2.3.2. in practice 

For the practical error computation, we perform a statistical dispersion of the observed 
points and assume a gaussian noise. The expression of the estimated noise for each 
observable is then the standard deviation of the measurements over an exposure: 


(Jx 


E" f (*-<x»2 


n 


2 

f 


( 2 . 11 ) 


3. Typical expected signal 

3.1. Simulating the companion 

We simulated a standard binary star model as in en. 13.1 l and applied the computation we 
currently use on the AMBER instrument to extract differential visibilities, differential 
phases and closure phase, computed the error bars assuming an average instrumental 
contrast of 50%, 1000 frames of integration, a detector noise of lie - per pixel and per 
frame and a total photon count of 4.2 x 10 7 , the same figures as in the following section. 
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CjkW 


1 + R(X)e~ 21 ^^ 
1 + R( A) 


(3-1) 


We used a modeled exoplanet spectrum from lBarman et all ( 2001 1 convolved with the 
resolution of the AMBER instrument in LR mode (R ss 35). We scaled this flux ratio 
to a maximum value of ss 10 -3 in order to estimate the signal amplitude of a 10 times 
brighter companion. 
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Figure 1. Simulations of an observation at VLTI with UT1-UT3-UT4 showing from left to 
right the differential visibility, the closure phase and the differential phase and from top to 
bottom separations ranging 1 to 3 mas. It shows a clear signal for all observables at 3mas but 
none for the closure phase at lmas whereas the visibility and phase still have detectable signal. 
Photon amount is 4.2 x 10 7 , detector noise lie - , number of frames 1000, number of pixels 32 
and average visibility 50%. These plots evidence the “super resolution” properties of differental 
phase and visibility relatively to closure phase. 


3.2. Results 

We tested the obtained signal for several Star/Companion separations and found that 
the signal detection at 3er would occur for a separation of only lmas for differential 
phase and visibility and 2.5mas for closure phase. This effects is due to the fact that if 
we perform a 1 st order taylor expansion of the phase, we get a linear dependence with 
ujk ■ ~ft whereas for visibility, we get a squared dependence and for the closure phase we 
get a cubic dependence. So for small separations, we have <^>diff > Vdiff > ip- 


4 . Best performances of AMBER 

We used the AMBER / VLTI instrument to observe the bright calibration star HD70060 
at low spectral resolution during the GTO run of 25 december 2004. We used the tech¬ 
nique explained below to compute the differential phases and differential visibilities and 
computed statistical error bars. These ones were compared to the theoretical ones as¬ 
suming a detector noise of 11 electrons and a null thermal noise. The figure 0 shows the 
resulting average differential phases for 5 successive exposures where we selected 50% 
of the best frames using the finge SNR as selection criterion. Each exposure represents 
about 20s of observation. They are separated by about 60 s. 

The standard deviation a in each exposure is the statistical dispersion of the differential 
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phase per spectral channel as described in eq. mu We performed weighted averages using 
the SNR on the fringe signal as weights, is the total number of collected photons per 
spectral channel and <T 0 is the accuracy expected from measured flux, detector noise and 
a supposed object visibility of 1 (which means the measured visibility is supposed to be 
only the instrumental one). 
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Figure 2. Differential Visibilities and Phases of the calibrator star HD70060 for 5 successive 
exposures (lines with letters) and the resulting average of the 5 exposures (thick line). 


4.1. Differential Visibility 

The differential visibility values vary only by 4 x 10 -3 with time within each exposure 
but the general slope of the curve changes dramatically between exposure, leading to 
a variation of about 0.05 radians rms over 5 minutes. This variation is dominated by 
the changes in the achromatic piston jitter due to seeing and vibrations fluctuations. 
In the current VLTI situation we have no tool to correct this effect, except including 
an estimation of the exposure jitter in the model fitting, with an impact on the SNR 
which cannot be estimated. So the differential visibility is currently not usable for very 
high accuracy applications. This situation will change dramatically when a fringe tracker 
is operational. Then we will be affected only by the residual piston jitter after fringe 
tracking correction. 


4.2. Differential Phase 

The rms variation of the measures within each exposure is of typically 1.8 milliradians (i.e. 
« 1/rarcsecond in colour-differential astrometry). This is about two times the expected 
rms from fundamental noise. Over the total 5 minutes we get 0.9 milliradians, again very 
close to twice the photon noise. This shows that the different measures seems statistically 
independent. An average of 1200 such exposures (15 hours) is needed to reach the 0.5 x 
10 -4 accuracy needed for the spectroscopy of r Bootis b. The brighter planet considered 
in m could be observed in 1 hour. 

With the improvement of the VLTI (less vibrations, improved overheads) we could 
expect to use almost all frames instead of only 50% of them with an average instrumental 
contrast improved by a factor 2. Then the r Boo observation would be achievable in a 
couple of hours. We also see a pattern as a function of lambda with a 10 ~ 2 radians rms 
over the K band. This pattern is stable over the 5 minutes of observations considered 
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here. That means that it can be eliminated by a fast calibration cycle. We think that 
it is a mixture of atmospheric dispersion and measurement effects. Measurement effects 
can be eliminated by beam commutation. Atmospheric dispersion will be eliminated in 
the closure phase and we plan to try to fit it in the differential phase. 

4.3. Closure Phase 

In this relatively poor quality early data, we have too little frames where the three 
fringe patterns are good enough. This explains why the closure phase is much more noisy 
(typically 10~ 2 radians rms) than the differential phases. We are therefore unable to say 
what part of the differential phase pattern is due to differential chromatic OPD. 


5. Conclusion 

The preliminary data reduction of bright sources observed in low spectral resolution 
with AMBER shows that the measured differential phases are accurate and stable enough 
to achieve the spectroscopy and angular separation of the most favorable Pegasi planets 
in a few 15 hours observations. This value should be reduced to 2 hours with the fore¬ 
seen simple improvements of the VLTI. The resulting spectra would be affected by an 
instrumental term and/or an atmospheric chromatic differential OPD term producing a 
smooth 10 -2 radians pattern over the K band. 

When the instrumental term will be eliminated by beam commutation, the remaining 
differential OPD might be possible to fit in the data reduction procedure. However, only 
a successful use of closure phase guarantees the elimination of the differential OPD. The 
current quality of the VLTI does not allow accurate closure phase measurements, but 
this should be improved soon, when the three fringe pattern are better stabilized. We 
remain very optimistic about the possibility to do spectroscopy of Pegasi planets with 
AMBER quite sson. 
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